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Abstract. Recently observational lower bounds on the strength of cosmic magnetic fields 
were reported, based on 7-ray flux from distant blazars. If inflation is responsible for the 
generation of such magnetic fields then the inflation energy scale is bounded from above 
as p\lf < 2.5 x 10~ 7 i\/pi x (i? o bs/10~ 15 G)~ 2 in a wide class of inflationary magnetogenesis 
models, where -B bs is the observed strength of cosmic magnetic fields. The tensor-to-scalar 
ratio is correspondingly constrained as r < 10~ 19 x (i? o b s /10 _15 G)~ 8 . Therefore, if the 
reported strength -B bs > 10~ 15 G is confirmed and if any signatures of gravitational waves 
from inflation are detected in the near future, then our result indicates some tensions between 
inflationary magnetogenesis and observations. 
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1 Introduction 

In 2010, the first detection of cosmic magnetic fields was reported [1] (see also [2-9]). Al- 
though High Energy Stereoscopic System (HESS) 7-ray telescopes observed TeV scale 7-rays 
from several blazars, Fermi space telescope did not observe GeV scale 7-rays from same 
blazars. Without cosmic magnetic fields, these two observations would contradict each other 
because some parts of TeV scale 7-rays traveling through the inter-galactic medium are con- 
verted into GeV scale 7-rays by electromagnetic cascade reaction; TeV scale 7-rays emit 
electron/positron pairs by scattering with extragalactic background lights and then created 
electrons and positrons emit GeV scale 7-rays by inverse Compton scattering with CMB 
photons until traveling C(l)Mpc typically [4]. In the presence of cosmic magnetic fields, on 
the other hand, they can bend the trajectory of charged particles and consequently decrease 
the flux of secondary GeV scale 7-rays. 1 From the observational lower limit on the bending 
angle the lower limit on the cosmic magnetic field strength was obtained. Several works were 
devoted to this subject [1-9]. Some of them took account of possible time variance of intrinsic 
blazar fluxes. The reported lower limit ranges O(10 _14 -10 _20 )G . 

The problem is what the origin of cosmic magnetic fields is. No astrophysical process or 
early universe phenomenology is known to explain sufficient amount of magnetogenesis [14]. 
As for the inflationary magnetogenesis, many scenarios were proposed, aiming to explain the 

1 Recently, refs. [10-12] contested the reliability of this mechanism, pointing out the possibility that the low 
flux of GeV scale 7-rays may not be due to cosmic magnetic field but may be due to a plasma beam instability. 
However, more recently, the Monte Carlo simulation in ref. [13] showed that the timescale of such instability 
is much longer than the timescale of electromagnetic cascade. Thus, according to them, the plasma beam 
instability cannot explain the lack of GeV scale 7-rays. 
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origin of magnetic fields of galaxies or galaxy clusters as well as cosmic magnetic fields [15- 
25]. The major obstacle in those models is the so-called "back reaction problem" [22, 26, 29]. 
Generation of magnetic fields during inflation inevitably increases the energy density of the 
electromagnetic field. If it becomes comparable with the energy density of inflaton then the 
dynamics of inflation is significantly altered. Consequently, the inflationary epoch may end or 
generation of magnetic fields is drastically suppressed. Demozzi, et al. [29] pointed out that 
in some specific models this problem is crucial and prevents generation of sufficient magnetic 
fields. 

In the present paper we conduct a model independent analysis of inflationary magneto- 
genesis. Specifically, we derive an upper limit on the inflation energy scale by assuming that 
all observed cosmic magnetic fields are generated during inflation. Our constraint depends 
on neither details of the model lagrangian, the behavior of photon mode functions nor the 
shape of magnetic field spectrum. If the strength of cosmic magnetic fields is stronger than 
10~ 15 G, the upper limit on the inflation energy density is 2.5 x 10~ Mp\. As a consequence, 
tensor-to-scalar ratio r is severely constraint as r < 10~ 19 . Therefore in this case, if all ob- 
served cosmic magnetic fields are generated during inflation, it is almost impossible to detect 
gravitational waves from inflation in near future observations. Conversely, if any signatures 
of inflationary gravitational waves are detected, then the cosmic magnetic fields should have 
another origin. 

The rest of the paper is organized as follows. In section 2, we introduce the observational 
lower limit on the strength of the magnetic field. In section 3, we discuss the assumptions 
which are needed to derive the main result. In section 4, the derivation of the upper limit 
on the inflation energy scale is presented. In section 5, we investigate the validity of the 
third assumption and explore the possibility to evade the constraint. Section 6 is devoted to 
a summary of this paper. In appendix, we derive the observational lower bound of cosmic 
magnetic fields in terms of the magnetic power spectrum. 



2 Observational Constraint on Magnetic Power Spectrum 

From the observations of blazars, current strength of cosmic magnetic fields is constrained 
[1-9] . While the constraints in those literatures are given in terms of the correlation length of 
magnetic fields, it is straightforward to rewrite it for the power spectrum of magnetic fields 
as 

^cfr(%ow) = J T F(kL)V B (r] now ,k) > B 2 ohs , (2.1) 
where VsiVnow, k) is the power spectrum of the magnetic field at the present, 



. . sin(z) . 
cos(z) V zbi(z) 



(2.2) 



and Si(z) denotes the sine integral function. (See appendix for derivation. For the introduc- 
tion of fediff) see the- explanation below.) Here, r\ is the conformal time, the subscript "now" 
denotes the present value, L = IMpc stands for the characteristic length scale for energy 
losses of charged particles due to inverse Compton scattering [4], -B b s is the observational 
lower limit on the strength of cosmic magnetic fields, and fcdiff — (1AU) -1 is the wave num- 
ber corresponding to the present cosmic diffusion length, i.e. the minimal size of a magnetic 
configuration which can survive diffusion during the Universe's lifetime [30]. Note that the 
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effect of the cosmic diffusion can be expressed as the cutoff in the power spectrum as 



V B (r) novn k)~0 (for k > fedjff). (2.3) 

in cases where magnetic fields have their origin in the primordial universe. In the formula 
(2.1), the cutoff at fc^iff is made more explicit as the domain of integration. 
For z > 0, F(z) satisfies 

< F{z) < 1, < zF(z) <a, a = Max[zF(z)] ~ 2.48. (2.4) 

Although -B bs still has a few orders of uncertainty, in the present paper we adopt the value 
reported by ref.[7] in which the latest data were analyzed. According to ref.[7], -B ODS — 10~ 15 G 
unless the time variance of intrinsic blazar fluxes is significant. 

The derivation of the formula (2.1) is given in appendix. Here, instead of showing the 
detailed derivation, we provide intuitive understanding of it. For this purpose, let us replace 
F(kL) by its asymptotic forms, 

F( Z ) ~ / 1 + °^ 2) (Z « 1} (2 V 



and drop 0(1) numerical factors to obtain the approximate formula as 



l / L dk 
k 



+ 



l/L k 



-j^PBiVnow, fe) 



(2.6) 



Let us now think of a Fourier mode of the magnetic field. For kL <C 1, the corresponding 
magnetic field can be treated as a homogeneous field, as far as the particle's trajectory (with 
the total length L) is concerned. Thus modes with kL <C 1 contribute to the bending angle 
as if they are homogeneous fields. This explains the first term in the right hand side of 
(2.6). On the other hand, for kL 3> 1, the direction of the corresponding magnetic field 
randomly changes N ~ kL times while the charged particle travels the total length L. If 
we were interested in the trajectory of the charged particle within one of short segments of 
the length ~ k~ 1 then the magnetic field could be treated as a homogeneous field. Actually, 
we are interested in the total bending angle due to N segments. Because of the randomness 
of the direction, the total bending angle from N segments adds up to only y/~N times the 
contribution from each segment. Therefore the contribution of modes with kL S> 1 to the 
variance of the bending angle should acquire the weight of order 1/N ~ l/(kL). This explains 
the second term in the right hand side of (2.6). 

3 Four Assumptions 

To derive the upper limit on the inflation energy scale, we need four assumptions. 
3.1 Assumption 1: the form of kinetic term 

First, we assume that the kinetic term of the photon field is of the form 

1 

4 



£kin = --firi)F^\ (3.1) 



where it is understood that the time-dependence of I(rj) is due to its dependence on homoge- 
neous, time-dependent fields present in the theory. Thus, £kin includes various interactions 
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between the photon field and other fields [16, 17, 22, 29]. This form of coupling does not 
have to break either gauge or local Lorentz symmetry. In general the photon field can have 
additional interactions £i nt : 

£-A = Adn + Ant- (3-2) 

However we let £; n t unspecified. Even so, under the four assumptions introduced in this 
section, we can derive the upper limit on the inflation energy scale in a model independent 
way. Note that when 1=1 and Ant = 0, the usual Maxwell theory is restored. 

This assumption on the form of the kinetic term is necessary to quantize the photon 
field and to define the kinetic energy density. The photon field can be separated into 
scalar and vector modes as 



AJr),x) = (Ao,Vi + diS) 



dM = 0, 



(3.3) 



where Aq and 5 are the scalar modes and Vi are the vector modes. Let us quantize the vector 
modes. After Fourier transformation with respect to the spatial coordinates, expansion by 
polarization vectors and mode expansion, we impose the standard commutation relation on 
the creation and annihilation operators. 



Vi(»7,x) 



(p) a m 
a k ' a k' 



p=i 



d 3 k 



ik-x^(p)/f> 



(2vr) 3 

(27r) 3 ^ 3 )(k-k')<5 pp ', 



t(p) u (p)*/ 



(3.4) 
(3.5) 



where Uk(rj) is the mode function of the photon vector mode, p (= 1,2) is the polarization 
label and e-^(k) is the polarization vector satisfying 



M P) (k) = (p = i,2), £^>(fc) e W(-fc) = * 



P =i 



(3.6) 



Then the canonical commutation relation for Vi requires the normalization condition of mode 
function u^rf) as 



4 P) 9^ { t 



(P)* Q JP) 



(p = l,2). 



(3.7) 



Now let us define the contribution of modes with k < &difT> i- e - those whose comoving 
length scale are longer than the cosmic diffusion length, to the kinetic energy density of the 
electromagnetic field as 2 



Pkin(&diff,f?) 



I 2 f kdia dk 

T 



[VE(v,k)+V B (v,k)} 



where we have defined power spectra of electric and magnetic fields as 
Ve{V, k) = — 2 b r , V B {V, k) - 



vr 2 a 4 (77) 



7r 2 a 4 (r?) 



(3i 



(3.9) 



2 It is understood that the domain of integration over k in eq.(3.8) is the same as in eq.(2.1). As a result, 
Pkin(fcdiff, rj) is not the sum of all modes that exist during inflation but the sum of modes whose scales are 
relevant to observed cosmic magnetic fields. 
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3.2 Assumption 2: Avoidance of strong coupling 

The second assumption is that 

I(rj) > 1 for VdiS < 7], (3.10) 

where r/difr is the conformal time when the comoving cosimic diffusion length exits the horizon, 
thus defined as fcdiff^diff = — 1- 

This assumption essentially states that the effective coupling constants of the photon 
field to other fields should be always smaller than present values. For example, let us consider 
the interaction between the photon and a charged fermion as 

Ant 3 -etirf^Ap. (3.11) 

In order to evaluate the effective coupling constant, we should canonically normalize the 
fields. Let us suppose that the fermion ip is already canonically normalized. The canonically 
normalized photon field is A^ = IA^. Then the interaction term is rewritten as 

Ant B -~^A%. (3.12) 

It is now clear that e/I is the effective coupling constant. Therefore if / <C 1, the effective 
coupling constant becomes large and the tree level analysis would be invalidated. In order to 
justify the tree level analysis, we need to assume that I is bounded from below by a positive 
constant Iq. For simplicity we set Jo to be the present value of /, i.e. Iq = 1. 

3.3 Assumption 3: Small back reaction 

The third assumption is that the kinetic energy density eq.(3.8) is smaller than that of inflaton 

3 

i 

Pkin^diff, V) < Pinf for r? diff < 77 < 77/, (3.13) 

where r/f is the conformal time at the end of inflation and hereafter we ignore the time- 
dependence of the inflaton energy density pj n f. 

This assumption is closely related to the condition for avoidance of the back reaction 
problem 

\Pkin(v) + Pint (V) I < Pinf for %iff <V<Vf- ( 3 - 14 ) 

Note that eq.(3.13) and eq.(3.14) are different. In general, the total energy density of the 
photon field includes not only the kinetic energy density p^ n but also the interaction energy 
density p- mt due to the additional interaction terms Ant- Also, pkin(^) in (3-14) should be 
understood as Pkm{oo,7]) and thus is in general larger than Pkin^diffi ? ?) i n (3.13). If the 
interaction energy density is non-negative {p- iri t > 0) then eq.(3.14) requires eq.(3.13). Even 
if the interaction energy is negative {p- mt < 0), unless the two contributions />kin arid pi nt 
cancel each other with a sufficiently good precision, eq.(3.14) generically requires eq.(3.13). 

3 We are interested in modes with k < fcdifi only. One can show that these modes do not significantly 
contribute to the vacuum enegy part of pkin- Validity of the effective field theory requires that the contribution 
of each mode uik to the vacuum energy must be smaller than the Planck mass scale Mpi. Hence, |p va c(fcdiff , f?) | < 
(Maxfc<fc d . ff \u)k |)(fcdiff /a) 3 < MpiH A <C p- in ( for Tjdiff < V < rjf. Therefore distinction between renormalized 
and unrenormalized expressions is irrelevant for (3.8). One can show that such distinction is unimportant also 
for (2.1). 
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Therefore the third assumption eq.(3.13) is mandatory unless negative pi Q t precisely cancels 

OUt positive /9kin- 

In section 5, we confirm the necessity of the third assumption in the case of gauge 
and local Lorentz invariant quadratic interactions and explore the possibility of the precise 
cancellation between p^ Q and pi n t. 

3.4 Assumption 4: Magnetogenesis during inflation 

The fourth assumption is that all observed magnetic fields are generated during inflation. In 
particular, the conformal symmetry of the photon field action is broken appreciably only in 
the inflationary era. Since the electric conductivity of the universe increases after the end of 
inflation [15] (see also (2.3)), we have 

5 c 2 fr(??now) < ajB 2 eS (rj f ), (3.15) 

where we have set a(r/ now ) = 1. 

By using eq.(3.10), (3.13) and the fact that B e g is smaller than the usual definition of 
magnetic field strength (0 < F(kL) < 1), we obtain 

B 2 cS (vdm) < 2pm. (3.16) 

Assuming the instantaneous reheating, we find the scale factor at the end of inflation is given 
by 

cA = (3.17) 

J Pint 

where p 7 ~ 5.7 x 10~ 125 Mpj ~ 5.2 x 10~ 12 G 2 is the present energy density of radiation. 
Eq.(3.15), (3.16), (3.17) and (2.1) lead to the following inequality. 

Iff^f < 10 " 2 x exph4(AAr - 35)1 (iFSb) 2 • (3 ' 18) 

where AN = ln(a f / a^is) ■ This inequality implies that 

4mB 2 cS ( Vdm ) « a}B 2 eS (vf), (3-19) 

and thus states that the magnetic fields have to be significantly amplified during inflation to 
explain the observational lower limit eq.(2.1). 



4 Upper Limit on Inflation Energy Scale 

With the four assumptions stated in the previous section, we are now ready to derive the 
upper limit on the inflation energy scale. The derivation is independent of details of infla- 
tionary magnetogenesis models, the behavior of photon mode functions or the spectrum of 
the electromagnetic fields. 

Independently from the specific functional form of Uk{rj), it can be shown that 

\ u k(Vf)\ 2 ~ K(%iff)| 2 = / dr? 2\u k (r])\ \u k {rf)\' 

2k\u k ( V )\ |4(77)| 

^k 2 \n k ( V )\ 2 + lu'^A , (4.1) 



< 



< 



Vf 
VdlB 

•If 

»7difl 



di] 
T 
drj 
~k 
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where we have used the general inequalities |-z(f?)|' < ^'(r/)! for a complex function z(rj) and 
2xy < x 2 + y 2 for real numbers x and y. Multiplying the both ends of eq.(4.1) by F(kL)k 4 /tt 2 
and integrating it over k from to fcdiff) we obtain 



cm fif f kdiB dk 

a}B 2 s (Vf) ~ 4 iS B 2 s (v d m) <j dr, a%) / - [V E ( V , k) + V B ( V , k)] , (4.2) 

where we have used the second inequality listed in (2.4). Using the second, third and fourth 
assumptions as well as eq.(3.19), we obtain 

B 2 s (r] now ) < ^Anf I d7]a 4 (ri) ~ ^" a)p inf (4.3) 



where H[ n f (~ const.) is the Hubble expansion rate during inflation. 

Note that 1/H m f and aj, which appear in the r.h.s. of (4.3), are decreasing functions of 
the inflation scale. Indeed, by substituting (3.17) for aj and using the Priedmann equation 
for -ffi n f, we can see that the r.h.s. of (4.3) is a decreasing function of the inflation scale. 
Hence, substituting eq.(2.1) and eq.(3.17) into eq.(4.3), we finally obtain the upper limit on 
the inflation energy scale, 

^ < ^ /4Mp1 ^ * 2 - 5 x 10 ~ 7Mpi x (iF^fe) " ■ (44) 

Note this upper limit can become even stronger if details of reheating is taken into considera- 
tion instead of eq.(3.17) [31]. Provided that the dominant energy density behaves like matter 
(oc a~ 3 ) during reheating, the right-hand side of eq.(4.4) is multiplied by an additional factor 
(/°reh/Pinf) 1//4 < 1; where /3 rc h is the energy density at the end of reheating era. 

Eq.(4.4) can be converted into the upper bound on the tensor-to-scalar ratio r under 
the slow-roll approximation, 

-19 I -Bobs 



Therefore, if all observed cosmic magnetic fields are generated during inflation, it is extremely 
difficult to detect any signatures of primordial gravitational waves, for example direct detec- 
tions or CMB B mode polarization. Conversely, if some observations reveal that r is larger 
than the upper bound (4.5), it implies that inflation cannot explain the origin of cosmic 
magnetic fields under the four assumptions. 

Now let us discuss the intuitive understanding of the reason why we obtain the upper 
limit on the inflation energy scale. Roughly speaking, c&Vb has to increase significantly 
during inflation for inflationary magnetogenesis (see eq.(3.19)). It is easy to show in the 
same way as eq.(4.1) that 

l^{a\ V )V B (v,k)) < ^(^(VEi^^) +V B (rj,k)) . (4.6) 

From the second and third assumption we know right-hand side of eq.(4.6) integrated by 
lnfe should be smaller than a pi n f ■ Thus essentially, time variation of ci^Vb is bounded from 
above by a 4 pj n f. Then we can rewrite a p; n f by using eq.(3.17) as 



a\r,)p- mi = Pl e- A \ N =I^L( ti -t). (4.7) 
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where N is e-folding number, t is cosmic time and tf denotes the end of inlation. Therefore 
a 4 pi n f is actually decreasing function of p in { during inflation. Since lower p m t is favored to 
relax the upper bound on time variation of o^Vb-, we obtain the upper bound on inflation 
energy scale. 

5 Additional Interaction Terms 

The action for the photon field consists of not only the kinetic term but also the ad- 
ditional interaction terms C[ n t- As already mentioned after (3.14), the third assumption 
eq.(3.13) is mandatory unless negative pint precisely cancels out positive pkin- Therefore 
whether such a precise cancellation is possible is a significant question. The answer we shall 
draw in the following discussion is that it is rather difficult to achieve such a cancellation. 

Here, it is perhaps worthwhile stressing that, as long as the four assumptions (including 
the third one) are satisfied, our main result eq.(4.4) holds even if pi nt and pkin precisely cancel 
out. 

5.1 Gauge and Lorentz invariant quadratic term 

In the quadratic level, the most general renormalizable interaction term which preserves 
gauge and local Lorentz symmetry is given by 

Ant = \f(n)e^F, u F pa + ^m\n)A,A^ (5.1) 

where e llvpa is the totally anti-symmetric tensor with e 0123 = g, f(rf) is a function of 

homogeneous scalars. The first term is called axial coupling term. The second term is the 
effective mass term of the photon induced by expectation values of charged scalars. It stems 
from the kinetic term of the charged scalars, and the positivity of the time kinetic term 
implies the positivity of the mass squared m 2 . This term spontaneously breaks the U(l) 
gauge symmetry, and the longitudinal mode of photon field becomes a physical degree of 
freedom. 

Actually, the axial coupling term does not contribute to the energy density of the photon 
field. Since the axial coupling term does not depend on the metric, its contribution to the 
energy momentum tensor is exactly zero. The effective mass term does contribute to pi n t but 
the contribution is always positive because of the positivity of the mass squared. Therefore 
the cancellation between pi nt and p^ n cannot occur. 

5.2 Model with negative interaction energy 

There is an existing model which gives a negative energy contribution from an additional 
interaction term. Turner and Widrow [15] proposed a model with non-minimal coupling, 
£; n t oc RA^A^, where R is the Ricchi scalar. This coupling can become an effective mass 
term of photon with negative mass squared. However this model has three critical problems. 
First, the longitudinal mode of photon becomes ghost [27-29]. Second, the negative energy 
contribution from £j nt exceeds p- m f and the back reaction spoils inflation when we require 
generated magnetic field is sufficient [29]. Third, this coupling explicitly breaks the gauge 
symmetry. 
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5.3 Energy conserving term 

From purely phenomenological viewpoints, let us investigate the additional interaction term 
of the form 

v^£ int = ^a 2 J 2 ( V )V? (J 2 > 0) (5.2) 

where J(rj) is a function of homogeneous scalar fields and Vi is the photon vector mode 
defined in eq.(3.3). This term is effective mass term of photon vector mode with negative 
mass squared. Note that this term has neither gauge invariance nor Lorentz invariance. It 
does not yield ghost field because it contains only vector modes by breaking Lorentz symmetry 
and we still assume that the kinetic term of photon is given by eq.(3.1). Although it may be 
hard to embed such a term in a viable elementary particle theory it is worth investigating it 
since we can find an interesting way to realize the cancellation between pi nt and /?kin- 
From eq.(5.2), the equation of motion is given by 

< + (k 2 - a 2 J 2 ) u k = 0. (5.3) 

Here we have assumed I{rf) = 1 for simplicity since otherwise the weak coupling effect due 
to / > 1 would make the interaction term irrelevant. At the same time we require the 
cancellation between p mt and p^ n for each mode, 

Kl 2 + (k 2 -a 2 J 2 )\u k \ 2 = 0. (5.4) 

It is easy to show that eq.(5.3), eq.(5.4) and eq.(3.7) imply that 

a 2 J 2 (r]) = const. (5-5) 

In other words, the coefficient of the quadratic term (5.2) should be constant. 

The reason why only the interaction term with constant coefficient leads the cancellation 
is simple. It is the energy conservation. If there is no explicit dependence on time in the 
action (for example, if the time-dependence due to the scale factor a(r]) is canceled by time- 
evolving scalars), then the energy of the system is conserved by virtue of Noether's theorem. 
In the case of eq.(5.2), if J(r/) cancels the time dependence of a(rj), the photon energy (with 
respect to the conformal time rf) is conserved. Note that the kinetic term of the photon 
field is originally free from a{rf). Therefore the energy density of photon does not increase 
even if the electromagnetic field strength increases. It is notable that, for this mechanism to 
work, the dynamics of the scalar fields included in J(rj) has to restore the time translation 
symmetry accidentally. 

The above analysis implies that the magnetogenesis from inflation whose energy is larger 
than the constraint of eq.(4.4) may not be impossible in principle. However, in practice it 
is not easy to realize a model which exploits the energy conserving mechanism because the 
accidental symmetry restoration by the scalar field dynamics can be easily spoiled by various 
effects such as the back reaction of the photon field. Therefore it is fair to say that all the 
four assumptions (including the third one) are likely to be mandatory in a rather broad class 
of models and the derived upper limit on the inflation energy scale is considerably general. 

6 Conclusion 

In this paper we have derived a universal upper limit on the inflation energy scale under 
the following four assumptions, (i) The kinetic term of the photon field is of the canonical 
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form up to a time-dependent overall factor, (ii) The effective coupling constants do not 
exceed present values and thus do not exhibit strong coupling, (iii) The kinetic energy of 
the photon field is always lower than the inflaton energy density during inflation, (iv) All 
observed cosmic magnetic fields are generated during inflation. 

The derived constraint is eq.(4.4), < 2.5 x 10" 7 Afpi x (S obs /10~ 15 G)" 2 . As a 
consequence, the tensor-to-scalar ratio r is bounded from above as eq.(4.5), r < 10~ 19 x 
(i? o bs/10 _15 G)~ 8 . We hardly expect that inflation is the origin of both cosmic magnetic 
fields and detectable gravitational waves if .Bobs > 10~ 15 G. Therefore the future detection of 
signatures of inflationary gravitational waves, if any, would imply tension between inflationary 
magnetogenesis and observations. 

Although our constraint is valid in fairly broad class of inflationary magnetogenesis 
scenarios, we have investigated the possibility to evade it. In order to evade the constraint, 
at least one of the assumptions should be violated. The third assumption can be violated 
only if the energy density due to additional interaction terms and the kinetic energy density 
precisely cancel out. We have considered a possible mechanism which exploits a energy 
conservation law to realize the cancellation. However, it seems a challenge to build a realistic 
model equipped with such a mechanism. 

So far, we have not succeeded to find a concrete model that produces relevant magnetic 
fields during inflation under the four assumptions we have made. However, this does not 
necessarily imply non-existence of such models. In this respect, it is probably worthwhile 
pointing out that our four assumptions do not exclude higher dimensional and/or nonlinear 
interaction terms. Our discussion can be applied also to models involving decays of other 
fields to the electromagnetic field. Those kinds of models have not been investigated in details 
and thus there is a large class of unexplored models. We feel that further discussion is beyond 
the scope of the present paper. 

Nonetheless, the result of the present paper is expected to be useful, providing a new 
judgment condition. Namely, if tensor-to-scalar ratio is detected in the future, any possibili- 
ties of magnetogenesis model within our assumptions will be excluded. Alternatively, if one 
can derive a lower limit on inflation energy scale by different arguments for a class of models 
then the upper limit obtained in this paper can be used to rule out the class of models. In 
this sense the upper bound we have found may be considered as an obstacle to inflationary 
magnetogenesis as well as an important guideline for model building. 
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A Derivation of the Constraint on Magnetic Power Spectrum 

In this appendix, we derive eq.(2.1). By Fermi and HESS observations, there is a lower 
limit on the bending angle of GeV scale cascade electrons and positrons in the inter-galactic 
medium. However, in the literatures the constraint on the cosmic magnetic field is given 
only in terms of the correlation length of magnetic fields [1-9] while theorists need the 
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constraint in terms of the magnetic power spectrum. In this appendix we shall generalize the 
constraint on the cosmic magnetic field to more general spectra. Such a generalization makes 
the connection between the cosmic magnetic power spectrum Vb and the bending angle 6. 
For simplicity, we neglect effects of special relativity in this appendix. 

Provided that a charged particle travels distance L in the background of a weak magnetic 
field B(r) from t\ till ti- Then the bending angle is given by 



e 



«(ti) -v(t 2 ) 



(A.l) 



where v(t) is the velocity vector of the particle. Note the absolute value of the velocity vector 
is constant. By using the equation of motion with Lorentz force, the difference of the velocity 
vectors is written as 



«(t 2 )-i>(*i) 



*2 



■/?? 



dt v(t) = — dt v(t) x B(t) = — dx x B(x) 



■m 



(A.2) 



where e and m are the charge and the mass of the particle, respectively, x(t) denotes the 
orbit of the particle and its initial value is set to x(t\) = 0. Then, we assume 9 is so small 
that the orbit can be approximated as a straight line, x(t) ~ x\(t)ei where e\ is the unit 
vector in the direction of the axis 1. By Fourier transforming B(x), we can perform the line 
integral 

d 3 k e iklL - 1 



dxi e\ x B[x\e\ 



(2vr) 3 iki 

By using these equations, we find that the variance of 6 is given by 



ei x B(k), 



(A.3) 



(o 2 



-)' 



2 fd 3 kd 3 k' i e 



\k\L 



1 



(27T) 



(5 i] -5 il 5 n )(B i (k)B j (k')), (A.4) 



Since the divergence of magnetic field vanishes (kiBi(k) = 0) and the cosmic magnetic fields 
are statistically isotropic and homogeneous, the square bracket in eq.(A.4) can be written as 



1 



(BiWBjik?)) = -(2tt) 3 8^ (fe + fe') 



k 2 J k 3 



2tt 2 

- 1 -^V B {k)+ie ijl k l H{k) 



(A.5) 



where Vsik) is the magnetic power spectrum and H(k) stands for the helicity component of 
magnetic fields [32]. By substituting eq.(A.5) into (A.4), we obtain 



(o 2 ) 



eL 
mv 



dk 
T 



V B (k) F(kL), 



3 2 



F(z) = -z 
v ' 2 



, . sin(z) 
cos(zj h zbi{z } 



l + 0(z 2 ) 

3tt 
Iz 



(z « 1) 



3n +0(z~ 2 ) (z>l) 



where Si(z) denotes the sine integral function. For z > 0, F(z) satisfies 

< F(z) < 1, < zF(z) <a, a = Max[zF(z)] ~ 2.48. 



(A.6) 
(A.7) 

(A.8) 



In order to find a proper definition of the effective strength of the magnetic field (in- 
cluding its normalization) for a given spectrum Vsik), as a fiducial configuration let us 
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consider a homogeneous magnetic field whose direction is perpendicular to the particle's tra- 
jectory. Denoting the strength of the fiducial magnetic field as B±, the bending angle is 
9 = L/Rl = (eB±L)/(mv), where i?L is the Larmor radius. On the other hand, for a sta- 
tistically isotropic spectrum, the variance of the magnetic field in three-dimensions is three 
halves of the variance of the magnetic field projected onto the two-dimensional subspace 
perpendicular to the particle's trajectory. Thus, it is natural to define the effective strength 
of the magnetic field as 



References 

[1] A. Neronov and I. Vovk, Evidence for strong extragalactic magnetic fields from Fermi 
observations of TeV blazars, Science 328, 73 (2010) [arXiv : 1006 . 3504]. 

[2] F. Tavecchio, G. Ghisellini, L. Foschini, G. Bonnoli, G. Ghirlanda and P. Coppi, The 

intergalactic magnetic field constrained by Fermi/LAT observations of the TeV blazar 1ES 
0229+200, Mon. Not. Roy. Astron. Soc. 406, L70 (2010) [arXiv : 1004 . 1329]. 

[3] C. D. Dormer, M. Cavadini, S. Razzaque, J. D. Finke and B. Lott, Time Delay of Cascade 
Radiation for Te V Blazars and the Measurement of the Intergalactic Magnetic Field, 
Astrophys. J. 733, L21 (2011) [arXiv : 1011 . 6660]. 

[4] H. Huan, T. Wcisgarber, T. Arlen and S. P. Wakoly, A New Model for Gamma-Ray Cascades 
in Extragalactic Magnetic Fields, Astrophys. J. 735, L28 (2011) [arXiv: 1106. 1218]. 

[5] K. Dolag, M. Kachclricss, S. Ostapchcnko and R. Tomas, Lower limit on the strength and 

filling factor of extragalactic magnetic fields, Astrophys. J. 727, L4 (2011) [arXiv: 1009. 1782]. 

[6] W. Essey, S. 'i. Ando and A. Kusenko, Determination of intergalactic magnetic fields from 
gamma ray data, Astropart. Phys. 35, 135 (2011) [arXiv : 1012 . 5313]. 

[7] A. M. Taylor, I. Vovk and A. Neronov, Extragalactic magnetic fields constraints from 
simultaneous GeV-TeV observations of blazars, Astron. Astrophys. 529, A144 (2011) 
[arXiv: 1101.0932]. 

[8] I. Vovk, A. M. Taylor, D. Scmikoz and A. Neronov, Fermi/LAT observations of 1ES 

0229+200: implications for extragalactic magnetic fields and background light, Astrophys. J. 
747, L14 (2012) [arXiv : 1112 . 2534]. 

[9] K. Takahashi, M. Mori, K. Ichiki and S. Inouc, Lower Bounds on Intergalactic Magnetic Fields 
from Simultaneously Observed GeV -TeV Light Curves of the Blazar Mrk 501, Astrophys. J. 
744, L7 (2012) [arXiv : 1103 . 3835]. 

[10] A. E. Broderick, P. Chang and C. Pfrommer, The Cosmological Impact of Luminous TeV 
Blazars I: Implications of Plasma Instabilities for the Intergalactic Magnetic Field and 
Extragalactic Gamma-Ray Background, Astrophys. J. 752, 22 (2012) [arXiv: 1106.5494]. 

[11] P. Chang, A. E. Broderick and C. Pfrommer, The Cosmological Impact of Luminous TeV 

Blazars II: Rewriting the Thermal History of the Intergalactic Medium, Astrophys. J. 752, 23 
(2012) [arXiv: 1106. 5504]. 

[12] C. Pfrommer, P. Chang and A. E. Broderick, The Cosmological Impact of Luminous TeV 

Blazars III: Implications for Galaxy Clusters and the Formation of Dwarf Galaxies, Astrophys. 
J. 752, 24 (2012) [arXiv: 1106 . 5505]. 




(A.9) 



Combining this with the formula (A. 6), we obtain 




(A.10) 



- 12 - 



[13] F. Miniati and A. Elyiv, Relaxation of Blazar Induced Pair Beams in Cosmic Voids: 

Measurement of Magnetic Field in Voids and Thermal History of the IGM, [arXiv : 1208 . 1761]. 

[14] A. Kandus, K. E. Kunzc and C. G. Tsagas, Primordial magnetogenesis, Phys. Rcpt. 505, 1 
(2011) [arXiv: 1007. 3891]. 

[15] M. S. Turner and L. M. Widrow, Inflation Produced, Large Scale Magnetic Fields, Phys. Rev. 
D 37, 2743 (1988). 

[16] B. Ratra, Cosmological 'seed' magnetic field from inflation, Astrophys. J. 391, LI (1992). 

[17] M. Gasperini, M. Giovannini and G. Veneziano, Primordial magnetic fields from string 
cosmology, Phys. Rev. Lett. 75, 3796 (1995) [hep-th/9504083]. 

[18] F. Finelli and A. Gruppuso, Resonant amplification of gauge fields in expanding universe, Phys. 
Lett. B 502, 216 (2001) [hep-th/0001231]. 

[19] A. -C. Davis, K. Dimopoulos, T. Prokopcc and O. Tornkvist, Primordial spectrum of gauge 
fields from inflation, Phys. Lett. B 501, 165 (2001) [Phys. Rev. Focus 10, STORY9 (2002)] 
[astro-ph/0007214]. 

[20] M. Giovannini and M. E. Shaposhnikov, Primordial magnetic fields from inflation???, 
[hep-th/0011105]. 

[21] B. A. Bassctt, G. Pollifronc, S. Tsujikawa and F. Vinicgra, Preheating as cosmic magnetic 
dynamo, Phys. Rev. D 63, 103515 (2001) [astro-ph/00 10628]. 

[22] K. Bamba and J. Yokoyama, Large scale magnetic fields from inflation in dilaton 
electromagnetism, Phys. Rev. D 69, 043507 (2004) [astro-ph/0310824]. 

[23] K. Bamba, C. Q. Geng and S. H. Ho, Large-scale magnetic fields from inflation due to 
Chern- Simons-like effective interaction, JCAP 0811, 013 (2008) [arXiv : 0806 . 1856]. 

[24] K. Bamba, N. Ohta and S. Tsujikawa, Generic estimates for magnetic fields generated during 
inflation including Dirac-Born-Infeld theories, Phys. Rev. D 78, 043524 (2008) 
[arXiv: 0805. 3862]. 

[25] R. Durrer, L. Hollcnstcin and R. K. Jain, Can slow roll inflation induce relevant helical 
magnetic fields?, JCAP 1103, 037 (2011) [arXiv : 1005 . 5322]. 

[26] S. Kanno, J. Soda and M. -a. Watanabe, Cosmological Magnetic Fields from Inflation and 
Backreaction, JCAP 0912, 009 (2009) [arXiv : 0908 . 3509]. 

[27] B. Himmetoglu, C. R. Contaldi and M. Peloso, Instability of anisotropic cosmological solutions 
supported by vector fields, Phys. Rev. Lett. 102 (2009) 111301 [arXiv : 0809 . 2779]. 

[28] B. Himmetoglu, C. R. Contaldi and M. Peloso, Ghost instabilities of cosmological models with 
vector fields nonminimally coupled to the curvature, Phys. Rev. D 80, 123530 (2009) 
[arXiv: 0909. 3524]. 

[29] V. Demozzi, V. Mukhanov and H. Rubinstein, Magnetic fields from inflation?, JCAP 0908, 
025 (2009) [arXiv: 0907. 1030]. 

[30] D. Grasso and H. R. Rubinstein, Magnetic fields in the early universe, Phys. Rept. 348, 163 
(2001) [astro-ph/0009061]. 

[31] V. Demozzi and C. Ringeval, Reheating constraints in inflationary magnetogenesis, JCAP 
1205, 009 (2012) [arXiv: 1202 . 3022]. 

[32] R. Durrer and C. Caprini, Primordial magnetic fields and causality, JCAP 0311, 010 (2003) 
[astro-ph/0305059] . 



- 13 - 



